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Abstract. We present results from an ab initio three-dimensional, multi-physics core 
collapse supernova simulation for the case of a 15 M e progenitor. Our simulation in¬ 
cludes multi-frequency neutrino transport with state-of-the-art neutrino interactions in 
the “ray-by-ray” approximation, and approximate general relativity. Our model exhibits 
a neutrino-driven explosion. The shock radius begins an outward trajectory at approxi¬ 
mately 275 ms after bounce, giving the first indication of a developing explosion in the 
model. The onset of this shock expansion is delayed relative to our two-dimensional 
counterpart model, which begins at approximately 200 ms after core bounce. At a 
time of 441 ms after bounce, the angle-averaged shock radius in our three-dimensional 
model has reached 751 km. Further quantitative analysis of the outcomes in this model 
must await further development of the post-bounce dynamics and a simulation that will 
extend well beyond 1 s after stellar core bounce, based on the results for the same 
progenitor in the context of our two-dimensional, counterpart model. This more com¬ 
plete analysis will determine whether or not the explosion is robust and whether or 
not observables such as the explosion energy, 56 Ni mass, etc. are in agreement with 
observations. Nonetheless, the onset of explosion in our ab initio three-dimensional 
multi-physics model with multifrequency neutrino transport and general relativity is 
encouraging. 


1. Introduction 


Numerous investigations over the past several decades clearly indicate the need to 

model core collapse supernovae in three sp atial dimensions wit h realist i c neut r ino trans- _ 

port an d general relativity [for a review, see lKotake et al.1 (l2QQ6l ): [jankal d2Q12l ): lMezzacappa et al. 
( 20151) 1, Realistic neutrino transport requires at a minimum a spectral, or multifre- 
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quency, treatment, where the lowest-order angular moments of the neutrino distribution 
function are evolved - e.g., the energy- and momentum-density per frequency. The 
models are therefore phase-space, and not purely spatial, models - i.e., at a minimum, 
carried out in the four-dimensional space of neutrino frequency and three spatial dimen¬ 
sions and, eventually, in the full, six-dimensional space of neutrino directions cosines, 
energy, and three spatial dimensions. Either of these approaches ("moments” or full 
Boltzmann) will require supercomputer architectures at the 100 PF - 1 EF scale, and 
perhaps beyond. Furthermore, the challenge to include a realistic treatment of neu¬ 
trino transport and general relativity adds significant complexity to both of these model 
components, of course, but especially to the neutrino transport dCardall & Mezzacappal 
12003k ICardall et all 12013bllah . Consequently, in an effort to advance to three-spatial 
dimensions with necessary physics - e.g., multifrequency neutrino transport and gen¬ 
eral relativity - one compromise has been developed and is used here: the ray-by-ray 
transport approximation dBuras et ahlfeOCGl) . In this approximation, a set of spherically 
symmetric transport problems are solved for each “ray.” A ray is an outgoing radial ray 
for each (6, p) pair - . Note, neutrino transport in this case is not restricted to pure radial 
transport. Rather, a complete spherical symmetric solve is performed. The net neutrino 
flux is however purely radial in this approximation. The more spherically symmetric 
the central source - in this case the proto-neutron star - the more valid the approxima¬ 
tion. A quantitative assessment of the quality of the approximation awaits an analysis 
using both ray-by-ray and non-ray-by-ray approaches in the context of an otherwise 
identical model. In the meantime, the ray-by-ray approach affords an ability to explore 
three-dimensions with the neutrino transport sophistication of past spherically symmet¬ 
ric models, at the expense of lateral neutrino transport. 

Recently, results from three-dimensional multi-physics models with detailed ray- 
by-ray neutrino transport and app roximate general relativity were reported by t he “Oak 
Ridge” and Max Planck groups ( Mezzacappa et al.l 120151: iMelson et ahll2015l) . Here, 
we focus on the results from our ongoing model. We report on the phases through the 
initiation of explosion. Additional analyses - e.g., a determination of the explosion 
energy - await the completion of 1-1.5 s of postbounce evolution. 


2. The CHIMERA Code 


CHIMERA is a parallel, multi-physics code built specifically for multidimensional 
simulation of core collapse supernovae. It is the chimeric combination of separate 
codes for hydrodynamics and gravity; neutrino transport and opacities; and a nuclear 
EoS and reaction network, coupled by a layer that oversees data management, paral¬ 
lelism, I/O, and control. The hydrodynamics are modeled using a dimensionally-split, 

Lagra ngian-Remap (PPM LR) scheme (IColella & Woodwardlll984h as impl emented in _ 

VH1 dHawlev et al.l2012h . Self-gravity is computed by multipole expansion dMuller & Steinmetz 
1 1995). In the GR case, the Newtonian monopole is replaced with a GR monopole 
(IMarek et al.ll2006l. Cas e A). Neutrino tra nsport is computed in the “ray-by-ray-plus” 

(RbR+) approximation ( Buras et~aD l2003h . Neutrinos are advected laterally (in the 6 
and <p directions) with the fluid and contribute to the lateral pressure gradient where 
p > 10 12 gcm -3 . The transport solver is an improved and upda ted version of the 
multi-group flux-limited diffusion transport solver of jBruennl (j 19851) enhanced for GR 
( Bruenn et al .112001 ). with an additional geometric flux limiter to prevent the over-rapid 
transition to free streaming of the standard flux-limiter. All 0(v/c ) observer correc- 
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tions have been included. CHIMERA solves for all three flavors of neutrinos and 
antineutrinos with four coupled species: v e , v e , v jJT = {v /( , v T }, v /IT = {vy,, v T ], with 
typically 20 energy groups each for ae = 4 - 250 MeV, where a is the lapse function 
and 6 is the comoving-frame group center neutrino energy. Our standard, modernized, 
neutrino-matter interactions include emission, absorption, and non-isoenergetic scatter- 
ingon free nucleons (iReddv et al .11 19981) . with weak mag netism corrections (Horowitz 


2002) : emission/absorption (electron capture) on nuclei (Langanke et al. 2003); isoen- 


ergetic scattering on nuclei, including ion-ion correlations; non-isoenergetic scattering 
on electrons and positrons; and pair emission from e + e~-annihil ation (lBruennlll985l) 
and nucleon-nucleon bremsstrahlung (lHannestad & RaffelH [l998 ). CHIMERA gener- 
ally utilizes the K = 220 MeV incompressibility version of the Lattimer & Swestv 
(Il99ll) EoS for p > 10 11 g cm' 3 and a modified version of the ICoopersteinl ( 1985 ) 
EoS for p < 10 11 g cm -3 , where nuclear statistical equilibrium (NSE) applies. Most 
CHIMERA simulations have used a 14 -speci es rr-nctwork ( a , 12 C- 60 Zn) for the non- 
NSE regions (XNet Hix & Thielemann 1999). To aid the transition between the net¬ 
work and NSE regimes, we have constructed a 17-species NSE solver to be used in 
place of the Cooperstein EoS for electron fractions Y e > 0.46. An extended version of 
the Cooperstein electron-photon EoS is used throughout. 

During evolution the radial zones are gradually and automatically repositioned 
during the remap step to track changes in the radial structure. To minimize restrictions 
on the time step from the Courant limit, the lateral hydrodynamics for a few inner zones 
are “frozen” during collapse, and after prompt convection fades the laterally frozen 
region expands to the inner 8-10 km. In the “frozen” region, lateral velocities are set 
to 0 and the lateral hydrodynamic sweep is skipped. The full radial hydrodynamics and 
neutrino transport are always computed to the center of the simulation for all rays. 

The simulation presented here utilizes 32,400 rays (solid angle elements) with 2° 
resolution in longitude and a resolution in latitude that varies from 8° at the pole to 
better than 0.7° at the equator, but is uniform in the cosine of the colatitude. Due to 
the Courant limit, the coordinate pole in standard spherical-polar coordinates creates 
a strong restriction on the time step size and therefore lengthens the total run time 
compared to a similar resolution two-dimensional simulation. Our constant cosine-of- 
colatitude grid seeks to minimize this impact without resorting to a grid that is coarse 
at all latitudes or implementing unevolved (frozen) regions near the pole. 


3. Results 


In Figured] we plot the angle-averaged shock radii as a function of time after bounce for 
our one-, two-, and three-dimensional models. All three models are initiated from the 


same 15 M e progenitor, taken from the most recent Woosley-Heger series (Woo sle v & Heger 
20071) , and include the complete physics set described in the previous section. After an 
initial leveling off of the shock radius at ~ 180 km between 75-100 ms after bounce, the 
outcomes begin to diverge, with the one-dimensional shock radius peeling off first and 
moving inward in a monotonic fashion until this particular simulation is halted. The 
two- and three-dimensional shock radii, however, maintain a value between 180 and 
220 km between 100 ms and ~ 190 ms after bounce, at which point the radii computed 
in the two- and three-dimensional models begin to diverge, with the shock radius in the 
two-dimensional case rising rapidly and the three-dimensional radius remaining at ~ 
200-210 km until it too begins to rise rapidly at ~ 280 ms after bounce. 
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Figure 1. Evolution of the shock trajectory from our ID model and the angle- 
averaged shock trajectories from our two- and three-dimensional models, all for 
the 15 M q case (iLentz et ak 201 5). The ID model does not develop an explo¬ 
sion, whereas an explosion is obtained in both our two-dimensional and our three- 
dimensional models. 


Entropy 



Time= 440.7 ms 



Z Axis (km) 


Figure 2. Snapshot of the equatorial cross section of the entropy in our on¬ 
going three-dimensional simulation for the 15 M 0 case at ~441 ms after bounce 
(ILentz et all 120151) . Red indicates high-entropy, expanding, rising material. 
Green/blue indicates cooler, denser material. Evident are significant (green) down 
flows fueling the neutrino luminosities. 
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Diagnositc energy 

3D v. 2D 



Figure 3. The diagnostic energy as a function of postbounce time for our ongoing 
three-dimensional model and its two-dimensional counterpart dLentz et al.ll2015l) . 


Figure 0 provides a snapshot of our ongoing three-dimensional model at 441 ms 
after bounce. Specifically, we show an equatorial slice through the postshock volume, 
of the stellar core entropy. The shock wave is outlined by the jump in entropy across it. 
Neutrino-driven convection is evident in the pattern of hot (red) rising plumes, which 
bring neutrino-heated material up to the shock, and cool (green) down flows, which re¬ 
place the fluid below. Clearly the distortion of the shock and the patterns of convection 
beneath it break axismmetry. Evidence for 1=1, “sloshing” and m = 1, “spiral” modes 
of the SASI will require a modal analysis, although the two-dimensional slice does not 
rule out either mode. Of particular note, at this instant of time one massive accretion 
funnel extends down to the proto-neutron star. Such funnels are a natural outcome of 
SASI-channeled postshock flows. Their importance lies in the fact they continue to fuel 
the neutrino luminosities on long time scales, of order 1 s after bounce. 

Figure [3] shows the evolution of the diagnostic energy in our three-dimensional 
model. Also plotted, for comparison, is the diagnostic energy for its counterpart two- 
dimensional model. The diagnostic energy is the sum of the total energy (gravitational 
plus internal plus kinetic) for all zones on our numerical grid for which this total is 
positive. The diagnostic energy does not take into consideration the work required to 
lift the stellar material above it, nor the energy gained in nuclear recombination as the 
initially outgoing material continues to move outward, potentially becoming unbound. 
With both of these included, we could define the explosion energy, but at this time after 
bounce, the explosion energy in our model is still not positive. Flowever, the diagnostic 
energy does in fact become positive at approximately 300 ms after bounce in our three- 
dimensional model. In the three-dimensional case, this is achieved at a time that is 
delayed relative to the time at which it is achieved in the two-dimensional case, by 
approximately 50 ms. 
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4. Conclusions 

Within the context of the approximations made, the results presented here indicate that 
the neutrino-driven explosions obtained in the context of two-dimensional core col¬ 
lapse supernova simulations also obtain here, in the context of a three-dimensional 
model. While the approximations made - most notably, the ray-by-ray neutrino trans¬ 
port approximation, as well as the general relativistic monopole correction to the New¬ 
tonian gravitational potential - must be assessed quantitatively in the context of future 
models that can accurately compare and contrast such approaches, we believe our use 
of detailed, spherically symmetric transport along each ray and a general relativistic 
monopole correction to the Newtonian potential in the context of the three-dimensional 
model presented marks a significant step forward. Clearly, an additional ~ 1 second of 
postbounce evolution is required in the model presented, and numerous models for dif¬ 
ferent progenitors, using different equations of state, etc. must be developed, as well. 
For these, we must have significant patience. Models such as the one presented here 
will run for a significant fraction of a year on either the DOE Oak Ridge Leadership 
Computing Facility or NSF’s Blue Waters platform. Thus, core collapse supernova 
modeling is entering a new phase, requiring the utilization of national-class resources 
for extended periods. As a result, the number of models that can be produced in a given 
year will be limited, and each model will require a significant time to complete. 
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